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Abstract 

The most general local Markovian stochastic model is investigated, for 
which it is known that the evolution equation is the Fokker-Planck equa- 
tion. Special cases are investigated where uncorrelated initial states re- 
main uncorrelated. Finally, stochastic one-dimensional fields with local 
interactions are studied that have kink-solutions. 

PACS numbers: 05.40.-a, 02.50.Ga 
Keywords: Fokker-Planck, stationary solutions, kinks 



1 mohamadi@alzahra. ac. ir 
2 mamwad@mailaps.org 



1 Introduction 



The Fokker-Planck equation was first used to describe the statistics of the Brow- 
nian motion of a particle in a fluid. This equation describes the time evolution 
of the probability density function of the position or the velocity of a Brownian 
particle. It is a differential equation, second order in the spatial variables and 
first order in time. There are also some studies on generalizations of the Fokker- 
Planck equation involving spatial derivatives of order higher than two [1-5]. 
Although these higher order equations are linear in the distribution function, 
these are sometimes referred to as nonlinear Fokker-Planck equations. There 
are also other generalizations of the Fokker-Planck equation [6,7]. In [8], it has 
been shown that any truncation of the generalized Fokker-Planck equation (the 
Kramers- Moyal expansion) is inconsistent, unless only terms of up to second 
order derivatives are kept. 

Reaction-diffusion systems is a well-studied area. People have studied reaction- 
diffusion systems using analytical techniques, approximation methods, and sim- 
ulations. A large fraction of exact results belong to low-dimensional (specially 
one-dimensional) systems, as solving low-dimensional systems should in princi- 
ple be easier. Despite their simplicity, these systems exhibit a rich and rather 
nontrivial dynamical and stationary behavior. Studies on the models far from 
equilibrium have shown that there is a remarkably rich variety of critical phe- 
nomena [9]. There are some exact results involving kinks in one-dimensional 
reaction-diffusion models, as well as simulations and numeric results [10], and 
also mean field results [11]. Formation of localized shocks in one-dimensional 
driven diffusive systems with spatially homogeneous creation and annihilation 
of particles has been studied in [12]. Recently, in [13], the families of models 
with travelling wave solutions on a finite lattice have been presented. Shocks 
have been studied at both the macroscopic and the microscopic levels and there 
arc efforts to understand how these macroscopic shocks originate from the mi- 
croscopic dynamics [11]. Hydrodynamic limits have also been investigated. In a 
recent article both stationary and dynamical single-kinks on a one-dimensional 
lattice have been investigated [14]. It was done for both an infinite lattice and 
a finite lattice with boundaries. Static and dynamical phase transitions of these 
models have also been studied. 

In such studies, one way is to investigate the system on a lattice and then 
take the limit of the results in the continuum. Another way is to write the 
equation in continuum from the beginning. In the latter case, one should write 
the master equation in continuum. This problem is addressed in this paper, and 
based on the results simple examples are studied. 

The scheme of the paper is the following. In section 2, the Fokker-Planck 
equation and its generalizions are briefly reviewed, where it is known that if the 
evolution equation of the probability density is local in space, then the evolution 
equation involves no more than second spatial derivatives [8]. By locality, it is 
meant that the time derivative of probability density is expressible in terms of 
only a finite number of spatial derivatives of that density. This means that 
the evolution equation for the most general local stochastic Markovian models 
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with is the Fokker-Planck equation. In section 3, this equation is generalized 
to the case of random fields. It is seen that there are special local interactions 
for which uncorrelated initial states remains uncorrelated. Finally, section 4 is 
devoted to some simple examples of stocahstic fields, admitting kink solutions. 

2 The Fokker-Planck equation 

To fix notation, let's briefly review the Fokker-Plank equation. Let r be a 
continuous (D-dimensional) random variable, and if)(r,t) be the corresponding 
probability density at the time t. One has 



V>(r,i) >0, 




(1) 



The evolution equation for this probability density is 





where U is the evolution operator with the following properties 



U(r,r';t,t')>0, 
d D r U(r,r';t,t') = 1, 



U{r, r';t,t) 



5{r-r'). 



(3) 



Differentiating ([2]) with respect to t, one arrives at 





or 





where H is defined through 





or 



H{v,v'-t) : 



dU(r, r'; 



t,t') 



(7) 



Using ©, one obtains 



H(r,r';t) > 0, 



r^r', 




(8) 
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Introducing the the basis vectors |r) and their dual basis covectors (r'|, 

(r'|r) =6(r-r'), 
Jd D r |r)(r| = 1, (9) 

one can write ([l} to ([5]) in a more compact form. To do so, one defines the 
vector and the operators H and U through 

ip(r) = (r\ip), 
H(t,t') = (r\H\r'}, 
U(r,r') = (r\U\r'}, 



and the covector (S\ through 



(S\ := / d^r <r|. 



One would have then 



>o, 

(Sty) =1, 

(r\U(t,t')\r>) >0, 
(S\U(t,t')=(S\, 
U(t,t) =1, 
d_ 

dt' 

au(t,t') 



\r/>(t)) =H(t) |V>(t)>, 



9i 



=H(t)U(t,t'), 
dU(t,t') 



dt 



t=t' 

(r\H(t)\r') >0, r^r', 
(5|fl"(t) =0. 



(10) 

(11) 

(12) 
(13) 
(14) 
(15) 
(16) 
(17) 

(18) 
(19) 

(20) 

(21) 
(22) 



These are in fact continuum analogs of a system with discrete states, and can 
be readily obtained from that by suitably changing summations to integrations. 

The evolution is said to be local (in space), if the time derivative of the 
density is expressible in terms of only a finite number of space derivatives of the 
density. We restrict ourselves to these cases. One then has 



ff-B-D**.-*, 



fc=0 



(23) 



fc=0 
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where f ll '" lk, s, and g ll "" lk, s are functions related to each other. This is equiv- 
alent to 



H(r,r';t) = £(-l) fc [d n ■ ■ -d lk 6(r - r')] (r'), 



fe=0 



: [^...^(r-r')] $p(r). 



(24) 



fc=0 



The fact that ^ is a probability density imposes constraints on its evolution 
equation. If tp(ro,to) is zero, then ip attains a minimum at ro and its first 
derivative vanishes at that point. As ip should remain nonnegative everywhere, 
one should have 



at 



> o. 



(25) 



The initial condition ijj(r,to) is arbitrary, except that it should satisfy ([T|). So 
the derivatives of order more than two of ip with respect to r at t = to and r = ro 
can attain arbitrary values. Suppose that of the derivatives (of order more than 
two) which appear in the right-hand side of (f2"3"|) , only one is nonvanishing and 
negative. If the absolute value of this derivative is large enough, then the right- 
hand side of (|2"3")) becomes negative, which contradicts (|2"5|) . So, one concludes 
that only derivatives of up to second order appear in the right-hand side of 
(|23|) . This is discussed in [8]. Moreover, at a point where ip vanishes, the 
matrix of its second derivative should be positive semi-definite, as at that point 
ip attains its minimum. The matrix g^) should be such that the time-derivative 
of ip be nonnegative at such a point. This means that g(2) should be positive 
semi-definite. 

So, the eveolution equation for ip is 

dip 
dt 
dtp 



■didi 



J(2) 



dt 



=5(0) i> ~ 9(i) dii> + 3 ( 2) didjip. 



(26) 
(27) 



Integrating (|2"6")1 on the whole space with ip(r, t) = 6(r — ro), and using (JTJ) , one 
deduces that /(o)(i"o) vanishes. As ro is arbitrary, this means that 

/(0) = 0. 



So (|26|) is simplified to 



dip 



+ didj 



(28) 



(29) 



Comparing this with (|27p , one obtains the relations between f(u\ 's and gn,) 's: 



= 



9(0) 



did 



/(i) + 2 9,9(2) ' 



J (2) -5(2) 



(30) 



4 



Introducing the operators R and P through 



it is seen that 



and 



R l \r) :=r*|r). 
Pi|r> ~di\r), 



\R\m =o, 

[R\Pj} =5), 
(SIP, = 0, 



(31) 



and the Hamiltonian for the system can be written as 
H=P i n i) {R) + P i P j fjUR), 



= - digfa + d t d 3 g^\ (R) + gfa (R) P t + gfa (R) P % P 

3 Random fields 



(32) 
(33) 

(34) 



Generalization of the arguments of the previous section to the case of fields 
is quite straight forward. A random field is a collection of random variables 
defined at each point of the space. One only needs to consider the following 
correspondences. 

k> He), 

R} -p(x), 

Pi -n(x), 

E- 

i 

dt- 



5q(x) 



(35) 
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where p is the random field operator, and one has 



p( x ) \q) : 


=q(k) 


Ie>, 


n(x) \ e ) ■- 


= Sg(^ 


[p(x),p(y)] : 


=0, 




[n(x),n(y)] = 


=0, 




[p(x),n( y )] : 


=<5(x- 


-y), 


(qW) ■ 


=S(g- 


-q'), 


(S\ :- 


= / Dg(g\, 


(5|n(x) : 


=0. 




It then follows that (|34|) is changed to 






H = J d n x n(x)/ (1) [x; p}+ J d n xd r 


l yn(x)n(y)/ (2) [x,y;p], 



(36) 



and /( 2 ) is positive semi-definite, that is, one has for arbitrary h and g, 

d n xd n yh(x)h(y)f {2) [x,y; Q }>0. (38) 



The evolution equation for ip (the probability distribution functional for p) is 
given analogues to ([26]) : 



d i =-f* n *sh {fwl * eMe]) 



+ / d ^ dny MxW (/(2)[X ' y ^ ] ^ ]) - ^ 

One may also obtain the evolution equation for the expectation value of the 
field ((p(z)) = (S\p(x\ip)). One has 

dp( ^ l I \ HI 
^b( z ),#l 



=/ (1) [z;p] + 2 /" d"xn(x)/ (2) [x,z;p], (40) 



from which 

|{p(z))^(/ (1) [z;p]). (41) 

It is seen that / 2 has no explicit contribution in the expectation value. It does, 
however, have implicit contribution through p. A special case is f% — 0. It can 
be shown then that if is a local functional of p, that is if [z; p] depends on 
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only p and a finite number of its derivatives in z, then if the system is initially 
uncorrelated in space, 



(p(x)p(y)) = (p(x))(p(y)}, x^y. 



(42) 



then it remains uncorrelated in space. For a system initially uncorrelated, and 
for x^y, one has initially 

d 

^(p( x )p(y)) =(/(i)[ x ;p]p(y)) + (p( x )/(i)[y;p]>, 

=(/( 1 )[x;p])( /5 (y)) + (p(x))(/ (1) [y;p]), 
(P(y)) + M x )> 



which shows that initially 
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[(p(x)p( y ))-(p(x))(p(y)}]=0. 



So the system remains uncorrelated. 
Applying the approximation 

(f(i){ x ,P\) « /(i)[ x , 
the evolution equation for the expectation value of p becomes 

d 



at 



(p(x))=/ (1) [x; (p)}. 



(43) 



(44) 



(45) 



(46) 



If /m is local, then the right-hand side of the above equation is the action of a 
differential operator on (p(x)}. 



4 Simple examples 

Consider a one-dimensional system for which is local, apply the approxima- 
tion (|45[) . and assume that 



(47) 



d 2 p 

fon x >P\ = + f Ip( x )}- 

Then the evolution equation for the expectation of p(x) is 

Of special interest are the time-dependent solutions to this equation: 

d 2 u(x) 



dx 2 



+ F[u(x)] = 0, 



(48) 



(49) 
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where 

u(x) := (p{x)). (50) 
Introducing the potential function V through 

^:=F(u), (51) 
one can once integrate (|49|) to obtain 

V(u) = const. (52) 
An obvious solution to this, or (149p . is 



2 \dx) 



u{x) = uq, (53) 

where uq is z zero of F. A kink solution is a solution which tends to constant 
values at infinity, but the value at +oo is different from the value at — oo: 

lim u(x) = u±, u + 7^ u_. (54) 

x — >±oo 

For such a solution to exist, one should have 

V(u-) =V(u+), 
dV 



du 



0, (55) 



V should be less than V(tt-) = V{u + ) if u is between u_ and u + . A simple 
example is 

V{u) = -a(u-ii_) 2 (M-u + ) 2 , (56) 

where u± and a are constants and a is positive. 

A multi-kink solution is a solution that tends to constant (different) values 
at plus and minus infinity, and is constant (or approximately constant) in some 
other intervals as well. Such solutions exist if there are some points Ui between 
u_ and u + , where the potential V becomes nearly equal to V(uJ) = V(u+). 
One then has a solution which tends to constants at infinity, and is nearly equal 
to Ui and nearly constant for a large region. An example is 

V(u) = -[a + P{u- u ) 2 } (u - u-) 2 [u - u+f, (57) 

where u±, a, and f3 are constants, a and (3 are positive, and uq is between u- 
and u+. If a is small, then there is an approximately double- kink solution. 

A special case is when V(uo) is exactly equal to V(u±). In this case, if V 
has second derivative in Uo, then the double-kink solution disappears and two 
distinct single-kink solutions emerge: one with between w_ and uq, the other 
between uq and u+. However, if the first derivative of V near uq tends to zero 
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like \u — u y with 7 < 1, then there is a solution with a flat part at u of a 
finite width. An example is 

V(u) =-0\u- u | 1+7 (u - u-) 2 (u - u+f, (58) 

where < 7 < 1. 

These are stationary solutions. It is easy to see that a stationary solution is 
readily converted to a travelling solution using a Galilean transformation. This 
transformation only changes the time-derivative Jj to + v • V, where v is the 
corresponding velocity parameter, and is in fact the velocity of the travelling 
solution. 
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